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Abstract The sum of factorials function, also known as the left factorial function, is de- 
fined as !n = 0! + 1!+--- + (n—1)!. These have been used by Smarandache 
and Kurepa to define the Smarandache-Kurepa Function (see reference [1], [2]). 
This paper presents an effective method for calculating !n, and implements the 
Smarandache-Kurepa function by using one new method. 


1. Introduction 


We define !n as 0! + 1!+---+(n—1)!. 

A simple PARI/GP program to calculate these values is below: 
ae 

sof f(n) = > i! 
i=0 

Then, 


for(i = 0, 10,print1(","soff(z))) gives the desired output; 
0,1,2,4, 10, 34, 154, 874, 5914, 46234, 409114, 
which is A003422 at OEIS [3]. 

2. A new method 


If we write out what the sum of factorials function is doing, we can write: 


1+ 
1+ 
1.2+ 
1.2.3+ 
1.2.3.4+ 
1.2.3.4.54+ 


and so on. 
If we now read down the columns, we see that this can be written as: 


LI [ip 2 3[L + 4[r be. 
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This is because we have an opening 1 from 0!. Then | is a factor of all the 
remaining factorials. However 1 is the only factor of 1 of the factorials, namely 
1!, so we have 


1+1[1+---. 

Having removed the 1!, 2 is now a factor is all remaining factorials, and is 

the final factor in 2!, hence 
1+1[1+2[1+--- 

and so on. 

!n requires inputs from 0! to (n — 1)!, and hence we are required to stop the 
nested recursion by n — 1. e.g. for !5, we have 

1+ 1[1+4 2[1 + 3[1 + 4{1]]]]. 


We can validate this: 


=141[1 + 2[1 + 36] 
=141[1 + 2[16]] 
= 1+ 1/33] 


= 34. 


3. Code for new method 


We can see how the new method decreases execution time, the original 
method presented performs O(k2) multiplications and O(k) additions. This 
method performs O(k) multiplications and O(k) additions. 

PARI/GP code for the routine is below: 

qsoff(n) = local(r); r = n; forstep@ = n —2,1,4-—l,re =t;r++);r 


4. Implementing the Smarandache-Kurepa function 


We need only consider primes, and the sk variable needs only range from 1 
to p— 1 (if !1 to !p are not divisible by p, then !(p + k) will never be as all new 
terms have p as a factor). 

For prime (p = 2, 500, for (sk = 1, p, if (qsoff(sk) 

This is obviously wasteful, we are calculating qsoff(sk) very repetitively. 
The code below stores the qsoffQ) values in a vector. 
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v=vector(500, z, qsoff(z)); forprime (p = 2500, for (sk = 1, p, if (v[sk] 
The following output is produced: 


2, 74, 0,0; 9,7, T5123 22516, a ames) 909, 04,42, -74,24,—5—,.20, 
86,2 507,= 183, =, = 64,,94,°70, 58: == A069, 1978 =, 14 = 958, 
167, —, 138, 80,59, —, —, —, —, 63, 142, 41, —, 110, 22, 286, 39, —, 84, —, —, 215, 80, 

14, 305, —, 188, 151, 53, 187, —, 180, —, —, —, —, 44, 32, 83, 92, —, 300, 16, —. 


5. Additional relations 


The basic pattern created in this paper also allows for the rapid calculation 
of other Smarandache-like functions based on the sum of factorials function. 

For example, we could define S'S F'(n) as the sum of squares factorial, e.g. 
SSF (10) = 0! + 1!+ 4! 4 9!, and the corresponding general expansion is 


UL O38 A 4 56.7-8.9/1 tee 


Or we can define the sum of factorials squared function as 
O!2+ 1124 2!24.---. 


In this case, the expansion is 
1+1[1+4[14+9/]1+---. 
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